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The gravitational ionization of a Keplerian binary system via normally incident periodic gravi- 
tational radiation of definite helicity is discussed. The periodic orbits of the planar tidal equation 
are investigated on the basis of degenerate continuation theory. The relevance of the Kolmogorov- 
Arnold-Moser theory to the question of gravitational ionization is elucidated, and it is conjectured 
that the process of ionization is closely related to the Arnold diffusion of the perturbed system. 
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1 I. INTRODUCTION 

Os\ n 

In a recent paperu, we considered the long-term nonlinear perturbations of Keplerian orbits by incident gravitational 
^"v waves of wavelengths much larger than the size of the system. In particular, we studied the periodic orbits of the 
perturbed system using the methods developed in Ref. [|. The existence of periodic orbits indicates the possibility of 
balance in the exchange of energy between the binary and the external radiation field. Thus gravitational ionization 
does not occur for such orbits. The issue of gravitational ionization is interesting as it involves the transport of energy 
by gravitational radiation and is analogous to the corresponding phenomenon that is well known in the electromagnetic 
context. 

I | The theoretical investigation of the interaction of a binary system with the gravitational radiation field reveals subtle 

■ phenomena that are further studied in this paper. In particular, the absorption of gravitational radiation energy by 
t-H ' the binary is not unidirectional in general. That is, the orbital energy of a binary immersed in a gravitational radiation 

field does not in general increase monotonically with time. On the other hand, the emission of gravitational radiation 
by the binary is expected to be accompanied by the monotonic decrease of the orbital energy of the system. In 
absorption, however, the incident wave can deposit energy into the orbit during one time interval and remove energy 
from the orbit during another time interval. A periodic orbit would result — even when the emission of the radiation 
by the binary is ignored — if after a certain time the net flow of energy between the incident wave and the binary is 

■ zero. 

^ Let us imagine, for the sake of simplicity, that gravitational radiation is incident on a Newtonian binary system 
I . consisting of a massive body of mass Mq at the origin of inertial coordinates and a particle of test rnaas to << Mo 
that revolves around it in the (x, y)-plane. The dynamical equation in this case is of the general formQu 



in 
o 



X 



dt 



•'" ^+efcy(t)x»'=0, (1) 



where k = Go{Mo + mo), e, < e << 1, is the perturbation parameter and e/Cy is the tidal matrix associated with 
the incident gravitational waves. Here JC is symmetric and traceless, and is related to the gravitational perturbation 
of the Minkowski space-time by 

Mt) = -i^(t,o), (2) 

where = i]^ + ex M „. Here we employ the transverse-traceless gauge for gravitational radiation, i.e. xo^ = 0, Xij 
is traceless and djXij = 0; moreover, Xij is a solution of the wave equation ^Xij — 0- Thus the gravitational radiation 
field can be expressed as a Fourier sum of monochromatic waves of frequency f2* and wave vector K*, c|K+| = fi*, 

Xij(t,x) = Re^ xV/( K *) ex P ' x - ito*t), (3) 
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where \ij is symmetric, traceless and XijK* = 0- The summation in (^) extends over all waves with much 
larger than the average orbital radius. Equation (Q) contains only the essential physics of the interaction of long- 
wavelength gravitational radiation with a Newtonian binary system; in fact, relativistic (i.e. post-Keplerian) effects 
in the binary are totally neglected. In particular, the emission of gravitational waves is ignored. The motivation for 
our treatment as well as its limitations is presented in detail in our recent workU 

The incident wave exchanges energy and angular momentum with the binary orbit but not linear momentum in 
the quadrupole approximation under consideration heraH. This is in exact analogy with the electromagnetic problem 
of the interaction of an electromagnetic wave with an atom in the dipole approximation. 

A simple linear perturbation treatment of (Q) has revealed the possibility of the existence of resonances at fi* = mu>, 
m = 1,2,3, • • •, where u> is the Keplerian frequency of the unperturbed elliptical orbit. Moreover, in this analysis 
secular terms appear that lead to the breakdown of the linear theory over timeQ. Thus linear perturbation theory is 
inappropriate for the investigation of periodic orbits of the perturbed system, since a periodic orbit is expected to 
persist forever beyond a certain point in time. 

In the first treatment of the nonlinear casetl, we considered a single monochromatic plane wave of frequency fi* that 
was normally incident on a Keplerian orbit of frequency u>. We found that in the generic case, certain orbits satisfying 
the resonance condition fi* = rruo.m = 1, 2, 3, • • •, could be continued to periodic orbits of the nonlinear system. The 
existence of periodic solutions of (|l|) demonstrates that ionization does not always occur; in fact, in a periodic orbit the 
energy exchange with the radiation field must be steady without any net flow. In addition, we found that for incident 
circularly polarized radiation of definite hclicity the rotation of the inertial coordinates by frequency fi*/2 rendered 
the dynamical equations autonomous. The invariancc of this autonomous system under time translation implies the 
existence of an energy integral in the rotating frame. In this case, the Kolmogorov-Arnold-Moser (KAM) theory 
implies that for sufficiently small e ionization can never occur in this system regardless of the magnitude of fi*/u;. 
To understand intuitively how this could come about, it should be pointed out that a binary system can gain or lose 
energy as it interacts with an incident gravitational wave. The situation in absorption is in contrast to the emission of 
gravitational waves by a binary. In the latter situation, the binary is expected to lose energy mjpnotonically; in fact, 
this is consistent with the observed rate of inward spiraling of the Hulse- Taylor binary pulsartffl. It follows that the 
reciprocity between emission and absorption of gravitational waves does not hold in general. This notion of reciprocity 
is valid in some other situations, however. For instance, a Keplerian ellipse of frequency to emits gravitational radiation 
of frequency mu>, m = 1, 2, 3, • • •, which corresponds to the resonance condition for absorption. 

Let us now consider a general periodic gravitational wave of period 2-7r/fi, t that is normally incident on the binary 
system. The existence of certain periodic solutions of the perturbed system may be expected on general grounds. It 
would therefore be more interesting to investigate the interaction of circularly polarized gravitational radiation with 
the binary system and to determine the stability of the resulting autonomous system under periodic perturbations. 
That is, a periodic wave may be expressed as a Fourier sum of components with frequencies nfi*, n = 1,2, 3, 
For a single component of definite hclicity, the transformation to the corresponding rotating frame would essentially 
remove the dependence of this perturbation upon time and for e below a certain limit the orbit would remain forever 
bounded even though there is a steady flow of incident gravitational radiation energy in the inertial frame. However, 
the time-dependence of the other Fourier components would not disappear in the rotating frame, and we would like 
to study the influence of these components on the ionization of the system. Though we develop methods that are 
applicable to a general periodic perturbation, we restrict our attention to a tractable problem for the sake of simplicity. 

In this paper, we consider a superposition of several harmonics in the perturbing function; clearly, the response 
of the system is not a superposition of the individual responses as a consequence of the intrinsic nonlinearity of the 
problem under consideration here. Specifically, we showed in Ref. [I] that for a normally incident circularly polarized 
monochromatic plane wave the motion is restricted to the (x, y)-plane and that a transformation to the uniformly 
rotating coordinate system in this plane with half the wave's frequency would result in an autonomous system for the 
equation of motion to which the KAM theory can be applied. It follows from the KAM theory that for sufficiently 
small e the motion is confined and ionization does not occur. We wish to explore the sensitivity of this interesting 
result to the particular form of the incident wave. Therefore, we consider here a principal right circularly polarized 
wave of frequency fi* = 2fi that is slightly modified by the presence of similar components of frequencies fi and 3fi 
as follows: 



Xu(t,0) = cos 2VLt + 25 
Xl2 (i,0) = sin 2fit + 25 



(a - (3) cos fit + i (a + (3) cos 3fit 
(a - ft) sin fit + i (a + 0) sin 3fit 



(4) 



X13 = X23 = X33 = 0. The other components of \ follow from the fact that x is a symmetric traceless matrix. Here 5, 
0<5<<l,isa new perturbation parameter that determines the relative strength of the extra secondary components 
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compared to the primary Fourier component of the normally incident radiation. Moreover, a and f3 are constant 
amplitudes of the order of unity, and the other numerical coefficients have been introduced for the sake of simplicity. 

The plan of this paper is as follows: In Section |3, we present the basic equations for a Hamiltonian description 
of the perturbed orbit in terms of Delaunay variables. Sections ^Il|-|v| are devoted to a development of degenerate 
continuation theory that is necessary for the identification of periodic orbits of the nonlinear problem via the methods 
and ideas that are originally due to Poincare. The existence of periodic orbits, described in Section |v|, demonstrates 
that a state of equilibrium can be established between the wave and the binary such that ionization does not occur; in 
fact, the net flow of energy vanishes in this case. To apply the KAM theory to our problem, it is best to transform (|l|) 
to a uniformly rotating frame as in Section VI . It follows from the description of the nonlinear system in this reference 
frame that Arnold diffusion is expected for S > 0, even for sufficiently small e. Numerical experiments described in 
Section VII tend to corroborate the conjecture that gravitational ionization is tantamount to Arnold diffusion in 
this system. For background material, this paper relies heavily on our previous detailed treatment of_the nonlinear 
problem for the case where the incident wave is essentially a simple monochromatic Fourier componcntEl; however, we 
have attempted to present sufficient detail here in order to render the present paper essentially self-contained. 



II. HAMILTONIAN DESCRIPTION IN DELAUNAY ELEMENTS 

Using (|l|) , (j^) and (Q) , we can write the associated Hamiltonian for this system as the sum of the Kepler Hamiltonian 
and the quadrupole perturbation given by ge/Cyrr'a;-'. The motion is taken to be in the (x, y)-plane, since the radiation 
is transverse and normally incident on the orbital plane; therefore, polar coordinates are convenient. Defining <f>(t) 
and ip(t) as 

4>(t) = n 2 {cos 2nt + -[(«- (3) cos fit + (a + 0) cos 3flt] } , 

iP(t) = n 2 {sm 2Q.t + - [(a - (3) sin M + {a + (3) sin 3fM]}, (5) 
the Hamiltonian in polar coordinates may be written as 

n = I fp' + 4) - - + er 2 [0(0 cos 20 + ?/>(*) sin 20]. (6) 
2 \ r z J r 

To express the Hamiltonian in a form particularly suitable for the analysis of orbital dynamics, we transform to 
the Delaunay elements (L,G,£,g). To this end, consider the bounded motion of the test particle according to the 
Hamiltonian (|6|). At each instant of time, the particle can be described as belonging to an osculating Keplerian ellipse; 
that is, the perturbed motion passes through an infinite sequence of osculating ellipses in the course of time. Each 
such ellipse is described by the unperturbed Hamiltonian 

H= 2\ P r + ^)--r' 

where we consider only bounded motions with E = H(p r ,pg, r,9) < 0. Let us now consider the canonical transforma- 
tion from (p r ,pg,r,9) to variables intrinsic to the ellipse. Thus we define action variables 

k 2 \ ^ 

L := [ —— , and G := pg, 



2E , 

which correspond to an osculating Keplerian ellipse with semimajor axis a, a = L 2 /k, and eccentricity e, 

L 2 

such that < e < 1. The equation of this ellipse is given by 

1 

r = a(l — e cos u), or r 



1 + e cos v 

where u is the eccentric anomaly and v is the true anomaly. The new canonical angle variables £ and g are then given 
by 
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I = u — e sin u, g = o — v. 

In the following, we exclude e = and focus attention instead on noncircular elliptical orbits. The resulting Hamil- 
tonian of the perturbed system in Delaunay variables is 

H = ~J^2 + 4C(L, G, i, g)<f>(t) + S(L, G, £, g)i>(t)], (7) 

where, in polar coordinates, C = r 2 cos 29 and S = r 2 sin 29. In Delaunay elements, C(L, G, £, g) and S(L, G, £, g) are 
given by 

C(L,G,£,g) = ^a 2 e 2 cos2 5 



-a 2 (A v (e) cos 2g cos vl — B v (e) sin 2g sin vl) , 



i/=l 

5 



S(L,G,e,g) = ^a 2 e 2 sm.2g 



oc 



+a 2 (Ay (e) sin 2g cos v£ + B v {e) cos 2g sin zW?) , (8) 



where 
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A v {e) = —(2ve(l - e 2 )J>e) - (2 - e 2 ) J„(i/e)), 



B v {e) = ~V^(eJiW - Kl - e 2 )J,(^ e )). (9) 

Here is the Bessel function of order v, and a prime indicates the derivative of the function with respect to 
its argument. The dynamical equations are derived in the usual way from the Hamiltonian (0). Moreover, these 
equations are given in Delaunay elements by 

• ( dC{L,G,£,g) ^ dS(L,G,e,9) , 

{ — qi — m + — oi — m 

■ ( dC(L,G,£,g) dS(L,G,e,g) 

{ Qg W) + ^ ^) 

dC(L,G,£,g) AM dS(L,G,£,9) , M 

dl m + di m 

dC(L,G,e,g) dS{L 7 G 7 £,g) , ,A n 
dG ^ + dG ^ ) ' (10) 

where uj is the Keplerian frequency of the binary given by to 2 = k/a 3 . 



III. BIFURCATION FUNCTION 



To establish the continuation (persistence) of periodic orbits of the Kepler problem to the system (10), we employ 
a method proposed in Ref. || and used in Ref. |l|; we only state the main ideas here and the reader is referred to these 
references for details. 



System ( |10[ ) has the abstract form 

u = F(u) + eh(u,t), (11) 

where u is a coordinate on a manifold M that consists of a cross product of Euclidean spaces and tori, and the 
function h is periodic with period 2ir/Cl in its second argument. We consider solutions t i— > u(t, £, e) of (nW with 
initial condition u(0,£,e) = £, £ € M, and define the mth order Poincare map by P m (£,, e) = u(27rm/f2, £, eJT Fixed 
points of this mth order Poincare map correspond to periodic solutions of (pd|). Consider the unperturbed periodic 
solutions of (O) with e = 0. These correspond to fixed points of the unperturbed mth order Poincare map defined 
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by p m (£) = 'P m (£, 0). Let us now suppose that there is a submanifold Z C M that consists of fixed points of p m , and 
( € 2. If there is a continuous curve e i— > «(e) in A/ such that re(0) = C and 7 ,m (K(e),e) = «(e), i.e. for each fixed 
value of e, «(e) is a fixed point of the mth order Poincare map, then the unperturbed periodic orbit is continuable. 
The continuation is established by the method of Lyapunov-Schmidt reduction to the Implicit Function Theorem 
and requires that every vector in TqM that is tangent to the submanifold Z be in the kernel of the infinitesimal 
displacement T>{Q = Dp m (Q — I. This is equivalent to the requirement that for each £ £ Z, the dimension of the 
kernel of the infinitesimal displacement at £ be equal to the dimension of the manifold Z. The manifold Z is called 
normally nondegenerate if it satisfies this condition. 

Let Z be a normally nondegenerate submanifold of M with dimension q; then, the range of the infinitesimal 
displacement at each point of Z has codimension q. For each ( £ Z, there is a vector complement <S(C) to the range 
of the infinitesimal displacement of dimension q. We denote the projection of T^M to 5(£) by s(£)- 

Let ( € Z and consider the curve in the manifold M defined by e ^ "P m (C, e); it passes through £ at e = and 
its tangent vector at e = is in TqM. This tangent vector can be identified with the partial derivative V™(C,0)- 
The bifurcation function B is defined to be the map from Z to the complement S of the range of the infinitesimal 
displacement 

s(0=S-(c)7T(C,o), 

so that in local coordinates B : R' -> R'. We define £ <E Z to be a simple zero of the bifurcation function if B(£) = 
and the derivative DB{Q is invertible. The following continuation theorem is proved in Ref. |2[ 

Theorem III.l Let Z denote a normally nondegenerate fixed point submanifold of M for the system fiT\). If ( G Z 

is a simple zero of the corresponding bifurcation function, then the unperturbed periodic orbit of with initial point 
£ is continuable. 



To apply Theorem III.l to the perturbed Kepler problem, we must compute the partial derivative P™(Ci 0) of the 
corresponding Poincare map. For the system (Eoh, the manifold M is the four dimensional Delaunay coordinate space 
and the Poincare map is defined as the strobe with period 2i:m/D, where f2 is the frequency of the perturbation. The 
partial derivative is obtained from the solution 1 1— > W(t) of the second variational initial value problem 

W = DF{u{t, C, 0))W + h(u(t, C, 0), t), W(0) = 0. 
In fact, we have W(t) = u e (t, (, 0) and 

7T(C,0) = W(2irm/n). 

We use Theorem [II. 1 to establish the existence of periodic orbits for the system (|l0|). Of course, this will establish 



the existence of periodic orbits for our model system ([[]), which describes the perturbation of a Keplerian binary 
system by a multi-frequency periodic gravitati onal w ave. We begin by identifying the normally nondegenerate fixed 



point submanifold of M mentioned in Theorem III.l . Recall that the frequency of the unperturbed periodic Keplerian 
orbit is u> — k 2 jl? . The three-dimensional manifold 

Z L := {(L,G,£,g) : mtu = nQ}, 

where m and n are relatively prime positive integers, is a normally nondegenerate submanifold of M (cf. Ref. [l]). 
Furthermore, the range of the infinitesimal displacement is complemented by the span of the vectors 



The bifurcation function associated with ( pp| ) is the projection of the partial derivative G, I, g, 0) on the 

manifold Z L onto the complement of the range of the infinitesimal displacement. To determine the bifurcation 
function, we solve the variational initial value problem 

L e = -^(L,G,£ + ojt, g)ct>(t) -^(L,G,£ + ut, fl )^(t), 

G e = -^(L,G,£ + wt, g)<j>(t) -^-(L,G,£ + ut, g)i>(t), 
og dg 

fin nc 

4 = -jjL £ + —{L,G,£ + ut, g)<t>(t) + —(L,G,£ + ut, g)^t), 
BC 8 S 
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with zero initial values; then, the solution evaluated at t = m(27r/f2) is projected to the complement of the range of 
the infinitesimal displacement. It follows from a detailed analysis that one can set the initial value of time equal to 
zero, as we have done, with no loss in generality. The bifurcation function is thus given by 

B(G,£,g) = (B L (G,£,g),B G (G,l,g),B*(G,£,g)), 

where 

r)T rff f)T 

B L {G,l,g):=-—, B G (G,£, g) := -— , B 9 (G,£,g) := — , (12) 

and 

f 2-nm/Q, 

[C(L, GJ + cut, g)<t>{t) + S(L, G,£ + ut, g)i>(t)] dt. 

To evaluate X, we use the resonance relation nil — mto to change the variable of integration from t to a = ttt/m+£/n, 
then observe that the integrand of X is periodic with period 2-k and substitute the Fourier series expansion for C and 
S. After performing these steps, we obtain the following expression for X in case n = 1 

X = irma 2 n{U 2m (e) cos (2g + 2m£) + -[(«- (3)U m {e) cos (2g + ml) 

+ (a + f3)U 3m (e) cos (2g + 3m£)}}, (13) 

U v (e) = A v (e) + B v {e). 



where 



If n = 2 and m is odd, then 

X = nma 2 QU m (e) cos (2g + m£); 
while for n — 3, if m is prime relative to 3, then 

X = ^irma 2 n6(a + (3)U m (e) cos (2g + ml). 

It turns out that X — for n > 3, as expected. 

We show in Ref. |l| that the bifurcation function does not have a simple zero for an incident monochromatic 
gravitational wave of definite helicity. The role of the secondary components of the wave is to resolve this issue for 
the bifurcation problem under consideration here. This is indeed the case for n = 1 as demonstrated in the next 
section. However, for n — 2 and n = 3 simple zeros do not exist, and the consideration of these cases would require 
the calculation of the solutions of higher order variational initial value problems. On the other hand, higher order 
perturbing functions of order e 2 , etc., are neglected in the formulation of equation ([j]), which is the starting point of 
our analysis. It follows from this remark that the treatment of the cases n = 2 and n = 3 is beyond the scope of this 
work. Thus, we assume n = 1 in the remainder of this section. 

Substituting ( p"3| ) into (|l2]), we obtain the following explicit form for the bifurcation function (n = 1) 

B L (G, i, g) = 2Trm 2 a 2 n{U 2m (e) sin (2g + 2m£) 

+-[(«- f3)U m {e) sin (2g + m£) 

+3(a + (3)U 3m (e) sin (2g + 3m£)}}, 
B G (G, I, g) = 2irma 2 n{U 2m (e) sin (2g + 2m£) 

+-[(«- 0)U m (e) sin (2g + m£) 

+ (a + /3)U 3m {e) sin (2g + 3ml)]}, 
G 



B"(G,£,g) = -Kma 2 n( — ){U' 2m {e) cos{2g + 2m£) 
+^[(a-0)U^(e)cos(2g + m£) 

+(a + /3)U^ m (e) cos (2g + 3m£)]}, (14) 



G 



where in the expression for B 9 we have used the fact that the eccentricity e and the Delaunay element G are related 



by G = ±LVT 



We therefore make a change of variable from G to e and observe that the zeros and their 



multiplicities for the bifurcation function B are identical to those of the function 

F(e,£,g) = (F L (e,£,g),F G (e,£,g),F 9 (e,£,g)), 



where 



F L (e, I, g) = U 2m {e) sin {2g + 2m£) 

+ -[{a - [3)U m {e) sin (2g + mi) 

+3(a + /3)U 3m (e) sin (2g + 3mi)], 
F G {e, £, g) = U 2m (e) sin {2g + 2ml) 

+ -[(a- (3)U m (e) sm(2g + m£) 

+ (a + p)U 3m (e) sin {2g + 3ml)], 
F»(e^ >fl ) = CC(e)cos(2 fl + 2mO 

+|[( a -/3)CC(e)cos(2 S + mO 

+(a + /9)^ m ( e ) cos ( 2 .9 + 3m£)]. 



(15) 



To apply Theorem [II. 1, we must determine the simple zeros of the bifurcation function ( |14| ) or, equivalently, (|1 
We will show that simple zeros exist by a perturbation argument which is presented in the next section. 



IV. ZEROS OF A DEGENERATE BIFURCATION FUNCTION 

The bifurcation function ( |l5| ) has the following abstract form 

AGu, e) = t(/*) +&p(M) = °> ^ r3 . (16) 
where the functions A : M 3 x I h> M 3 and r, ip : R 3 i— > R 3 are given in components by 

/ t i(m) \ / ¥>iO) \ 

r(/i) = r 2 (/x) , ip{fi) = ^ 2 (m) , 
\ t s(m) / V ^ 3 ^) / 

and where the first two components of r are equal. In fact, we define r = t\ = t%. 

Remark IV. 1 The fact that the bifurcation function has the abstract form ( fl^ j is not accidental. Every normally 
incident gravitational plane wave whose Fourier representation has as its dominant term a purely circularly polarized 
wave will result in a bifurcation function of the form ( f7^ ) at resonance. 

We will determine the zero set for A in case e is sufficiently small. 

Lemma IV. 2 Suppose A is defined by (fjfi)). If r) € R 3 is such that T{rj) = while the vectors gradr(?7) and gradr3(7y) 
are linearly independent, then there is a curve s > T(s) in R 3 such that T(0) = i], r'(0) ^ 0, and A(T(s),0) = 0. If 
such a curve exists and s = is a simple zero of the real-valued function s i— » ip 2 (r(s)) — tpi(T{s)), then there is a 
curve e~ i— » T(e) m R 3 swc/i £/ia£ T(0) = rj and A(T(e), e) = 0. Moreover, for each sufficiently small t ^ 0, the point 
T(e) €E R 3 is a simple zero of the function fi t— * A(/x, e). 

Proof. The linearly independent vectors gradr(?7) and gradr3(?7) span a plane Q in R 3 . Let S denote one of the two 
possible rotation operators in space that preserves Q and rotates each vector in Q through ir/2 radians. We note that 
Sgradr^) and SgradTa (77) are linearly independent vectors in Q. Also, we let v denote a vector in space such that 
the set {Sgra,dT(r]), Sgia,dT 3 (r]),v} is linearly independent. This requirement is satisfied, for instance, if we choose v 
to be equal to the cross product of Sgr&dT(r]) and iSgradT3(?7). 
The function w : R 3 i— > R 3 defined by 

((ii,(4i,H3) !-> A*i'S'g r adr3(77) + n 2 Sgra,df(r]) + /j 3 v + r) 
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is invertible. We use it to define 9 : M 3 — > M 2 given by 

9(Ml,M2,M3) = C^M/Jl,/^,/^)), T 3 (w(m,fl 2 , M3))- 

It follows that O(0) = 0. Moreover, the derivative of the transformation 

^ M2,0)),t 3 (w;(/xi, /i 2 ,0)) 

at (121,112) = (0,0) is given by the matrix 

grad?(?7) • Sgradr 3 (?7) 

gradT 3 (r7) • Sgradr^) 



(17) 



that represents the partial derivative of O with respect to its first two arguments at the origin. 

Using the linear independence of gradr(?7) and gradr3 (77) , it is easy to see that the diagonal elements of the matrix 
( p^7| ) are both nonzero. Thus, by the Implicit Function Theorem, there is a unique curve s 1— > <r(s) in the (/xi, ^2)-plane 
such that cr(0) = (0,0) and 6(ct(s),s) ee 0. The curve T(s) := w(a(s),s) is such that T(0) = rj and A(r(s),0) = 0. 
Moreover, 



T'(O) = Dw(0) ( CT '[ 0) 



Since Dw(0) is invertible and the vector (cr'(0), 1) ^ 0, we have that T'(0) 7^ 0. This proves the first assertion of the 
theorem. 

Under the assumption that s = is a simple zero of the function s 1— > (p2(T(s)) — <px(T(s)), we will use the Lyapunov- 
Schmidt reduction procedure to show that the corresponding point rj = T(0) in K 3 is continuable as a zero of A. As 
r'(0) ^ 0, the image Z of T is locally a one dimensional submanifold of R 3 . Moreover, Z is normally nondegenerate 
in the sense that at each point z G Z sufficiently close to 77, the kernel K of Dt(z) is exactly the one dimensional 
tangent space of Z. In fact, since this tangent space is clearly in K, it suffices to show that K is one dimensional. 
The derivative of r is expressed in the standard vector notation by 

gradf (z) 
Dt(z) — I gradr(z) 
gradr 3 (J) 

Since the vectors gradr(z) and gradr 3 (;j) are linearly independent at z = 77, they will, by continuity, remain linearly 
independent in an open neighborhood of r\. Thus, their span is two dimensional at each point of this neighborhood. 
In particular, for z in this neighborhood, it is clear (by matrix multiplication) that gradr(z) and gradr 3 (z) are not 
both in K. This proves that K is one dimensional. It follows from the same Lyapunov-Schmidt reduction procedure 



that is used to prove Theorem [II. 1 (see Ref. g) that if TL(z) denotes a projection to the complement S of the range 
of Dt(z), then simple zeros of the map from Z to K given by z 1— > II(z)Aj(z, 0) are continuable. Thus, if 77 is such a 
simple zero, there is a curve e'l— » T(e) such that T(0) = 77 and A(T(e),e) = as required in the lemma. 

To construct such a projection into S, it can be shown that, in fact, the range of Dt{z) is spanned by the 
transpositions of the vectors (1,1,0) and (0,0,1). To see this, it is sufficient to note that -Dt(z) Sgradr (z) is the 
transpose of a scalar multiple of (0,0,1) while £>T(z)S , gradT 3 (z) is a scalar multiple of the transpose of (1,1,0). 
The transpose of the vector (0, 1,0) clearly spans a complement to the range that we will denote by S. Thus, if 
(pi,p2,pa) G M 3 , the projection II is easily computed and is given by II(px, P21 Pa) = P2 — Pi- Note that II does not 
depend on the base point zeZ. 

Using the projection II and the definition of A, the Lyapunov-Schmidt reduced function 2 i— > II(z)Ag(z, 0) is given 
by z 1— > <fi2(z) — fi{z) f° r % G Z. Thus, if, in the coordinates of Z, s — is a simple zero of s ip 2 (T(s)) — ipi(T(s)), 
then 77 = r(0) is a continuable zero of A given by a curve T, as required in the statement of the lemma. 

It remains to show that T(e), for sufficiently small et^ 0, is in fact a simple zero of the function /1 1— > A(/j,e). To 
this end, it suffices to show that the matrix DA(T(e),e) does not have a zero eigenvalue. But, under the hypotheses 
of the lemma, the matrix DA(t7,0) has a one dimensional kernel, namely the tangent space of Z. This means that 
DA(r], 0) has exactly one zero eigenvalue. By the continuity of eigenvalues of matrices, there is a smooth family 
of eigenvalues A(e) such that A(0) = and a corresponding smooth family of eigenvectors V(e) such that V(0) is a 
nonzero vector tangent to Z with 
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DA(T(e),e)V(e) = X(e)V(e). 

It suffices to show that if e is positive and sufficiently small, then A(e) ^ 0. By continuity, the remaining two 
eigenvalues of the matrix will be nonzero as well. The desired result follows as soon as we show that the derivative 
A'(0) ^ 0. For this we have 

D 2 A( V , 0)(T'(0), V(0)) + DAg(r), 0)1/(0) + DA(r h 0)V"(0) = X'(0)V(0), 

where D denotes differentiation with respect to the space variable (i€l 3 . After projection by IT into the complement 
of the range of DA(r], 0), 

nz) 2 A(?7,o)(r'(o),y(o)) + nx>A e -(r/, o)v(o) = A'(o)nv(o). 

We claim that UD 2 A(rj, 0)(T'(0), V(0)) = 0. Once this claim is proved, the fact that A'(0) ^ would follow 
provided HDAi(r], 0)V(0) ^ 0. The key point to note is that the projection operator, II, is independent of the base 
point, i.e. the range of the projection operator is spanned by a constant basis vector independent of the base point. 
Thus since A e -(T(s),0) = <p(T(s)), we have 



= nxv(r7)r'(o) = uDA^(rj, o)r'(o). 



s=0 



Since Z is one dimensional, T'(0) is just a scalar multiple of V(Q) and therefore if HD(p(r])T' (0) ^ 0, then 
UDAi(ri, 0)V(0) 0. By the definition of the components of ip and of the projection II, we have 



d j / \ \ 

-n„(r (s)) 



s=0 aS 



s=0 



Thus, if s = is a simple zero of s <— > (p2(T(s)) — tpi(T(s)), then 

ODAg(»7,0)V(0) ^ 0, 

as required. 

To verify the claim, note that since II projects to the complement of the range of DA(rj,0) and since II does not 
depend on the base point on Z, we have HDA(T(-d), 0)T'(#) = for the real variable $ in the common domain of T 
and T. This implies that 

0= A[n£>A(r(tf),0)T'(tf)] 

aV i?=0 

= nL> 2 A(77, 0)(r'(0), T'(0)) + n£>A(?7, 0)T"(0) 
= mD 2 A(r,,0)(V(0),r>(0)), 

where cV(0) = r'(0). Since c ^ 0, the lemma is proved. 



V. PERIODIC ORBITS 



We now use Lemma IV. 2 to prove that some of the zeros of the bifurcation function 
that once this result is established 



[III 



15| ) are simple. We note 
that the corresponding bounded orbits of the 



it will follow from Theorem 
Keplerian two-body system are continuable under perturbation by periodic gravitational waves 
If 5 = 0, then the zeros of (|l5|) are the union of the following one dimensional sets 

Z+ = {(e, £, g) : 2g + 2m£ = 0, U' 2m {e) = 0}, 
Zr = {(e,£,g) : 2g + 2ml = tt, U' 2m {e) = 0}, 

Z+ = {{e 1 i 1 g):2g + 2m£=^, U 2m (e) = 0}, 



Z 2 ={{e,e,g) : 2g + 2ml = 



2' 
3tt 

T' 



U 2m (e) = 0}. 



(18) 







We will show that the zeros in the set s Zf~ and Z 2 continue to simple zeros of the bifurcation function for sufficiently 
small (5 7^ 0. To conform with Lemma IV. 2 , we use ( |l5| ) to identify the components of A as they appear in the lemma 
as follows: 

r(e, I, g) = U 2m (e) sin (2g + 2ml), 
r 3 (e, I, g) = U' 2m {e) cos (2g + 2ml), 

(f!(e,l,g) = - [(a- f3)U m {e) sin {2g + ml) 
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-3(a + (3)U 3m {e) sin (2g + 3ml)} 



(p 2 (e,l,g) = - [(a- 0)U m (e) sin (2g + ml) 



2 



+(a + (3)U 3m (e) sin (2g + 3ml)} , 



<p 3 {e, l,g) = - [(a - /3)f4(e) cos (2g + ml) 



2 



+ (a + l3)U^ m (e) cos (2g + 3ml)} 



Also, we note that 5 in (|15| ) plays the role of e in the lemma. 

There are four cases to consider corresponding to the zero sets ( |l8| ) of the unperturbed bifurcation function (Jl4|) . 
We will consider the zero set for illustrative purposes, as the computational procedure is identical in all four cases. 
The set Z± is just a line in K 3 , so here the curve T in the lemma can be taken to be a parametrization of this line 
starting at an appropriate point. By a direct calculation, the function ip 2 — <p± on Z^ is given by 

- [(a - f3)U m (e) + (a + (3)U 3 m (e)] sin g. 

The zeros of this function along Z^ are clearly simple provided the coefficient of sin g does not vanish at the value of 
the eccentricity determined by membership in Z^ . The main result of this section is the following theorem. 

Theorem V.l Consider the system and suppose that [L, G,l, g) is on an unperturbed (m : n) resonant periodic 
solution with n = 1, that is, L 3 — mAr/fi. If (e,l,g) is in Z^ or Z^ , respectively, (e,l,g) is in Z 2 or Z 2 , where e 
is the eccentricity of the corresponding Keplerian ellipse (e 2 = 1 — G 2 / L 2 ) and if 

(a - (3)U m (e) + {a + f3)U 3m (e) £ 0, 

respectively, 

(a - (3)U m (e) - (a + p)U 3m (e) £ 0, 
then the Keplerian ellipse continues to a periodic orbit under the perturbation. 

It remains to show that the sets {Zf }, i = 1,2, are not empty. This fact will follow as soon as we show that both of 
the functions Ui m and U 2rn have zeros on the interval < e < 1. 

Recall that U2 m (e) — A 2m (e) + B 2m (e) and note that both A 2rn and B 2m have removable singularities at e = 0. 
Moreover, both have Taylor series at e = with leading terms given by 

2m -2 

In particular, if m = 1, then lim e ^ + [/2(e) = 2 and for m > 1, the limit is zero, but U^mifi) > for sufficiently small 
eccentricity. Also, linie^x- ^2m = — J2m(2m)/m 2 . By a standard property of the Bessel functions, J v (u) > 0; hence, 
Uimil) < 0. This proves U2 m has at least one zero on the interval < e < 1. Numerical calculations suggest that 
this zero is unique. 

A simple argument can be given to prove the existence of a zero in the case m = 1 for the function U 2m . In this 
case, the series expansion of U 2 at e = is — lOe + 0(e 2 ). Thus, U 2 (e) < for a small but positive eccentricity. The 
limit as e — > 1~ is the same as 

e->l~ vl — e z 

By standard properties of the Bessel functions, Ji(fS) — and J2 > on the interval (0, f 2 -J, where f 2 x is the first 
zero of J 2 . Hence, J' 2 > on (0,j' 2 x ). But, we also have 2 < j' 2 1 . Thus, J' 2 {2) > and U 2 (e) — * 00 as e — > 1 _ . This 
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proves that U' 2 {e) has at least one zero on the interval < e < 1; moreover, numerical computations suggest that this 
zero is unique. 

For m > 2, we will outline a proof that shows the function U' 2m has at least two zeros on the interval < e < 1. 
An asymptotic analysis shows that the function U' 2m is positive near the end points of the interval. The proof is 
completed by showing that U 2m has a negative value within the interval. To this end, note that the function U 2m has 
the form 

u 2m( e ) = li(jn,e)J 2m (2me) + 72 (m, e) J' 2m {2me), 

where 71 and 72 are computed using (ft). Moreover, by standard properties of the Bessel functions, both of the 
functions e 1— > J2 m (2me) and e 1— ► J 2m \2me) are positive for < e < 1. Let us define e := (1 — l/(4m 2 )) 1/ ' 2 , and 
observe that < e < 1. An simple computation shows that 

7i(m, e) < 0, 7a(m,e) < 0. 

Hence, we have U 2m (e) < 0, as desired. Numerical computations suggest that the function U 2m has exactly two zeros 
on the interval < e < 1. 

It follows from these results that the resonant interaction of the incident multi-frequency gravitational wave of 
definite hclicity with a Keplerian binary can result in orbits that are periodic with period 2ir/Q, where ti — mu>, 
m = 1,2,3, . . .. It is conceivable that other periodic orbits may exist; however, our method can identify only those 
periodic orbits that are continuations of resonant Keplerian orbits of the unperturbed system. 



VI. ROTATING FRAME 

In a manner similar to that used by Hill in his treatment of the lunar theoryETH, we can view the dynamical system 
described by (|^) in Cartesian coordinates rotating at half the principal frequency of the incident gravitational wave 
(O* = 20). These coordinates — that we again represent by (x,y) — rotate with frequency \ti± = ti with respect to 
inertial coordinates; therefore, the equations of motion in these coordinates are given by 

fl - 2ti^ -fi 2 x + ^- + 2eti 2 [(1 + 6a cos tit) x 
at 2 dt r d 



(80sintit)y] = 0, 

fy + 2fl— - n 2 y + % + 2eQ 2 [(60 sin fii) x 
dt 2 dt 



-{l + 8acosQ£)y\ =0. (19) 

Let X :— x — Qy and Y := y + Qx be the canonical momenta conjugate to x and y, respectively. Then ( |l9| ) is 
equivalent to a Hamiltonian system with Hamiltonian 

U R = \{X 2 + Y 2 ) + n(yX -xY)-- 
2 r 

+eQ 2 [(l + Sacosnt)(x 2 - y 2 ) + 2(60smQt)xy\. (20) 

By identifying the momenta in polar rotating coordinates as p r = (xX + yY)/r and pg — xY — yX, ( p0[ ) in polar 
coordinates is given by 



H R = - [pi + p J - - - Qp e + en 2 r 2 [cos26 

+6(a cos 29 cos fit + sin 29 sin tit)], (21) 
which in the corresponding Delaunay elements becomes 

k 2 

U R = -— 2 - tiG + eti 2 {C{L, G, i, g) + 6[aC(L, G, i, g) cos tit 

+0S(L,G,l,g)smtit]}. (22) 



The Hamiltonian system given by (|22|) has periodic orbits with period 2ir/ti, where ti = mu>, m = 1,2,- • •; this 
assertion can be demonstrated using results that have already been obtained in this paper. In fact, using the resonance 
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assumption, periodic orbits exist in the inertial frame with period 2ir/Q as shown in Section An orbit in the inertial 
frame, t i— » (xi(t), yz(t)), is represented in the rotating frame by 

XB.{t) =xi(t) cos fit + yi(t) sin fit, 
Unit) = —Xi(t) sin 17i + yi (t) cos 17i. 

It follows from these relations that if xi(t) and yi(t) are periodic with period 27r/17, then so are Xa(t) and yn(t). 
Therefore, the orbits that are periodic in the inertial frame with period 2ir/Q are observed to be periodic with period 
27r/rZ in the rotating frame. It is also possible to arrive at this conclusion by direct application of the methods 
developed in the previous sections. 



To understand the physical structure of Hamiltonian (22), imagine right circularly polarized gravitational radiation 
of frequency 217 that is normally incident on the orbital plane. It follows from our previous workEJ, as well as equation 
©, that in a reference frame rotating with frequency 17 the wave stands still. That is, observers at rest in the 
rotating frame do not perceive the variability associated with a wave so that in the absence of secondaries (6 = 0) 



the dynamical system (19) is autonomous. That an observer — by merely rotating about the propagation axis of the 
circularly polarized wave — could make the wave stand completely still would be a remarkable physical effect and 
deserves further discussion. 

It is a fundamental consequence of Lorentz invariance that all basic radiation fields travel with speed c with respect 
to all inertial observers. This may be illustrated by an example: Let an inertial observer move with speed vq along the 
propagation axis of a monochromatic plane gravitational wave of frequency 17*. The frequency and the wave vector 
of the radiation as perceived by the moving observer are smaller than those measured by static inertial observers by 
a common Doppler factor of 

- vp 
+ v 

Mathematically, as v — > c this ratio goes to zero and hence the frequency and wave vector of the radiation vanish 
so that the wave might appear to stand still. This limit is not physically allowed, however. No observer can move at 
the speed of light, although — theoretically-p=pne can get arbitrarily close. Therefore, the wave can never stand still 
for an inertial observer. It has been shownli!3 that according to the standard Einstein theory this is not the case for 
accelerated observers, i.e. an accelerated observer can indeed stand still with respect to a gravitational wave. The 
autonomous nature of the system (|l^) for 6 = provides an interesting illustration of this fact. That is, consider an 
observer at the center of a system of coordinates rotating with frequency 17. The observer does not move, but the 
fact that it refers its observations to the axes that rotate with frequency 17 with respect to the inertial axes makes it a 
noninertial observer. Radiation of frequency 17* is incident in the inertial frame along the axis of rotation. According 
to the noninertial observer, the frequency of the gravitational wave is 17^ = 17* =F 20, where the upper sign refers to 
right circularly polarized (RCP) gravitational radiation and the lower sign refers to left circularly polarized (LCP) 
gravitational radiation. The first (second) case has helicity +2 (—2), so that 17^ = 17* — h • f2, where h is the helicity 
of the gravitational radiation field; this is an example of the general phenomenon of helicity-rotation coupling. Now if 
17 = ^17* for RCP gravitational waves or 17 = — ^17* for LCP gravitational waves, we find that 17^ vanishes according 
to the noninertial observer (as well as any other observer at rest in the rotating frame anywhere along the z-axis) and 
the radiation field stands still. It has been shown in Ref. [I] that for a monochromatic gravitational wave with definite 
helicity and sufficiently small amplitude this observation concerning a rotating observer leads to the conclusion that 
a Keplerian system in the presence of this radiation can never ionize. In principle, this absence of ionization could be 
considered an observable consequence of the physical possibility that a gravitational wave could stand completely still. 

Let us now consider the possibility of ionization of the binary system as a function of the parameter 6. The 
transformation to the rotating frame leaves the orbital radius unchanged; therefore, the ionization problem can be 
discussed equally well in the rotating frame. In fact, the problem becomes simpler since the principal component of 
the incident radiation field loses its time - dependence as in (|2^). Thus if 8 = 0, the KAM theorem implies that for 
sufficiently small e the perturbed trajectory is bounded since it is trapped between two-dimensional invariant tori in 
the three-dimensional energy surfaces. When <5 7^ 0, however, the two secondary components in the inertial frame 
both reduce to a perturbation of frequency 17 in the rotating frame; that is, the three Fourier components of the 
radiation field in the inertial frame are RCP waves with frequencies 17* = 17, 217 and 317, while the tidal matrix in the 
rotating frame has frequencies given by 17^ = 17* — 217, i.e. —17, and 17 for the three copipanents, respectively. The 
perturbation of frequency 17 in the rotating frame is expected to lead to Arnold diffusioiiE3II3 and hence ionization of 
the system. However, it has not been possible thus far to prove ionization for the system (|2^); therefore, we resort to 
numerical work in the following section. 
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VII. NUMERICAL EXPERIMENTS 



In this section we illustrate — by means of numerical experiments — the conjecture that gravitational ionization and 
Arnold diffusion are closely related. The interpretation of the numerical results is simplified if we take the viewpoint 
of inertial observers and consider the perturbed motion given by Hamiltonian (||). 

We have performed several numerical experiments to test the diffusion and ionization properties of the dynamical 
system that is represented by the Hamiltonian (^3j). The physical meaning of these numerical experiments is essentially 
the same as in our previous paper (cf. Figure 2 in Ref. 0]) : Let us choose two scales for the measurement of time and 
length that are arbitrary except that they are connected here by our choice of k = 1 . In these otherwise unspecified 
units, we have chosen an unperturbed ellipse of semimajor axis a — 4/3 and eccentricity e = 1/2 such that g = —n/2 
and the Keplerian frequency is u = 3\/3/8. The ellipse is perturbed by the presence of periodic gravitational radiation 
and the orbital radius of the osculating ellipse is then plotted versus time in these units in Figures 1-4. These figures 
represent the results of our numerical experiments in which we have set the parameters of the external perturbation 
as follows: 



and have changed S in the range < 5 < 1 and f2 in the range 1 < Q < 2. In each run, we set the initial conditions 
to be 



which correspond to the unperturbed ellipse described above. After integration over each time interval corresponding 
to one cycle of the perturbation, i.e. 2ir/fl, the corresponding elapsed time t and orbital radius are plotted. Clearly, 
for each value of t (abscissa) there is only one value of orbital radius r (ordinate); however, this is not discernible 
in some of the figures due to the way in which the plots have been prepared. The KAM theorem is illustrated in 
Figure 1, where the value of r appears to oscillate between ~ a(l — e) = 2/3 and ~ a(l + e) = 2, thus indicating the 
complete absence of ionization. The possibility of ionization of the system is illustrated in Figures 2 and 3 for 5 > 0, 
where the initial ellipse is near resonance in the top panel, on resonance in the middle panel and off resonance in the 
bottom panel. More precisely, the middle panel in either Figure 2 or Figure 3 corresponds to an exact third-order 
resonance since fl/to = 3, while the top panel illustrates the response of the ellipse to the external perturbation near 
resonance Cl/ui 3.08; the bottom panel is off resonance with fl/ui « 1.54. Additional calculations extending the 
integration time for the system depicted in the top panel of Figure 3 have been performed. These results suggest that 
a bursting behavior occurs in which the near resonance condition shown at the right end of this panel is followed by 
chaotic motion similar to that shown in the middle of this panel that in turn is followed by a period of near resonance. 
This recurrence of chaotic and near resonance behavior appears to continue for the extended interval of time studied 
as illustrated in Figure 4. Thus, this behavioij4s consistent with a type of chaotic behavior, called intermittent chaos, 
that has been studied for dissipative systemai3. If the behavior suggested by these simulations is indeed present in 
the Hamiltonian system (|^), then our result would be an example of Hamiltonian intermittency. 

Imagine, for the sake of concreteness, a binary system consisting of an artificial satellite in an eccentric orbit 
about the Earth. Let the scales of length and time be i?o and To, respectively; then, i?g = fcTg . Thus, if we take 
i?o = 10 9 cm for the problem under consideration, it turns out that To 1.6 x 10 3 sec. The gravitational wave in our 
numerical experiments would then have a frequency of the order of fl « 10 -3 rad sec -1 , corresponding approximately 
to 1.5 x 10~ 4 Hz as well as to a wavelength of 2 x 10 14 cm, and an amplitude of the order of e = 10~ 3 . Gravitational 
waves have not yet been directly observed; however, in a realistic situation the amplitude of the wave would be 
expected to be of the order of 10" 20 . 

Our numerical experiments are consistent with the expected behavior for a 2 ^-degree of freedom Hamiltonian 
system. In fact, by introducing a fictitious action variable (cf. §6 in Ref. [I]), our system is equivalent to an autonomous 
Hamiltonian system with three degrees of freedom. Each orbit of this new system is constrained to an energy manifold. 
However, whereas the two dimensional KAM tori (if they exist for our choice of parameter values) separate each three 
dimensional energy manifold for the corresponding two-degree of freedom Hamiltonian system that we obtain with 
6 = 0, the three dimensional KAM tori, that may exist for the three-degree of freedom Hamiltonian system that we 
obtain with S ^ 0, do not separate space within the five dimensional energy manifolds. Of course, for sufficiently small 
choices of 6, there are orbits of the three-degree of freedom Hamiltonian system that remain bounded for all time; for 
example, the periodic orbits of Section |y| and the orbits confined to KAM tori. While the totality of bounded orbits 
in an energy manifold may be a set of positive measure, we expect that every open set of the five dimensional energy 
manifold contains an initial condition for a trajectory that will diffuse throughout the energy manifold. In fact, we 
expect this behavior for all 5 > 0. On the other hand, as 8 decreases toward zero, the time required to leave the 
vicinity of a KAM torus is expected-,to grow at an exponential rate. Figures 1-4 illustrate dynamical behavior that 



e = .001, a = 



2.5, (3 = 2 



(p r ,pe,r,0) = (-5,1,1,0), 
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VIII. APPENDIX 



The purpose of this appendix is to point out that the main results of this paper still hold for a more general incident 
gravitational wave than that considered in (0) . 

In the general case of incident gravitational radiation on a binary system, the motion of the system away from 
the initial orbital plane (i.e. along the ^-direction) also needs to be taken into account. It is necessary to mention 
here that initial conditions are generally ignored in our theoretical approach, which relies on the properties of the 
perturbed system once transients have died away and a "steady state" situation has been established. In this paper, 
we have limited our considerations to normally incident waves; in fact, the transversality of gravitational radiation 
makes it possible to set z = for normally incident waves. We wish to note here that essentially the same results can 
be obtained for a more general incident radiation field. 

Consider, for instance, the superposition of a left circularly polarized (LCP) wave of frequency fl traveling along 
the z-axis with a linearly polarized wave of frequency fl traveling along the x-axis. Specifically, let Xij be of the form 



cos fit — sin fit 
X(t, 0) = ( sin fit 3 cos fit 

o 0-4 cos fit 



(23) 



up to a constant factor. It follows from (g) that the corresponding JC is of the form 

jc = ^fi 2 x (t,o), 

so that the equation of motion ([!]) along the z-direction is given by 



(24) 



dt 2 



^ - 2efl 2 z cos fit 



0. 



(25) 



This equation is satisfied by z — 0, which is consistent with the fact that the orbit is always in the (x, y)-plane. 
Therefore, to the incident radiation field in (0J one could add radiation fields of the type given by (23). Moreover, 
upon transformation to the rotating frame 



1 



cos fit sin fit 



JC' = R^JCR = -fl 2 I sinflt 3 cos fit 
1 




(26) 



where 



cos fit — sin fit 
R = | sin fit cos fit 
1 



(27) 



is the rotation matrix used in Section VII. It is remarkable that the tidal matrix JC' in the rotating frame has the 
same frequency fl as in the incrtial frame; in fact, JC' can be obtained from IC simply by letting fl — > — fl. 
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FIG. 1. Orbital radius versus time plots for the dynamical system given by Hamiltonian (|6j) with initial conditions 
(p r ,pe,r,0) equal to (.5,1,1,0) and parameter values e = .001, 8 = 0, a = 2.5, /3 = 2, k = 1, and, from top to bottom, 

n = 2, Q. = 9V3/8, n = 1. 



FIG. 2. Orbital radius versus time plots for the dynamical system given by the Hamiltonian (|6|) with initial conditions 
(p r ,pe,r,8) equal to (.5, 1, 1,0) and parameter values e = .001, <5 = 0.5, a = 2.5, j3 = 2, k = 1, and, from top to bottom, Q = 2, 
n = 9^3/8, Q = 1. 



FIG. 3. Orbital radius versus time plots for the dynamical system given by the Hamiltonian (|6|) with initial conditions 
(p r ,pe,r,0) equal to (.5, 1, 1,0) and parameter values e = .001, 5 = 1, a = 2.5, /3 = 2, k = 1, and, from top to bottom, Q — 2, 

n = 9^3/8, n = i. 



FIG. 4. Orbital radius versus time plots for the dynamical system given by the Hamiltonian (|6|) with initial conditions 
(p r ,pe,r, 6) equal to (.5, 1, 1, 0) and parameter values e = .001, 5 — 0.5, a = 2.5, /3 = 2, k = 1, fi = 2. This is an extended form 
of the plot presented in the top panel of Figure []. 
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